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Abstract 

Given a finite set X of distinct points, Marinari-Mora's 'Axis of Evil 
Theorem' states that a combinatorial algorithm and interpolation enable 
to find a 'linear' factorization for a lexicographical minimal Groebner 
basis (J(I(X)) of the zerodimensional radical ideal /(X). 
In this work we provide such algorithm, showing that it ends in a finite 
number of steps and that it actually provides the correct result. 
The 'Axis of Evil' algorithm takes as input the monomial basis of the 
■ initial ideal T(I(X)) but its starting point is the (finite) Groebner escalier 

N (obtained via Cerlienco-Mureddu correspondence) so we will also define 
the 'potential expansion' 's algorithm, a combinatorical algorithm which 
computes the minimal basis from a finite Groebner escalier. 



00 

U 
< 



; 1 Introduction. 



> 
in 

ON 
OO 



X 

5-H 



Keywords: Groebner basis , Combinatorial algorithm, Interpolation. 



Marinari-Mora in [10], [9], [TT] gave a deep description of the structure of a 
zero-dimensional ideal / described by giving its Macaulay basis B(I) ([S]); in 
' particular they enhanced the description of the Grobner basis of an ideal in 

P\) . K[X, Y] given by Lazard in [5] proving that in a restricted case which includes 

the radical one, for each monomial r := X^ 1 ■ ■ ■ X£™ belonging to the mini- 
mal basis G(I) of the initial ideal of /, it is possible to produce linear factors 
7m5r := Xi — f(X\, . . . , 1 < m < n, 1 < 5 < d m such that the polynomi- 

als f T := Hm=i Il5=i ImSr form a minimal lexicographical Groebner basis of I; 
each such factors were obtained by producing an appropriate decomposition of 
the given Macaulay basis B(I) = Um=i U<5=i SmsM and interpolating over the 
monomial set obtained applying Cerlineco-Mureddu Algorithm over the set of 
functionals S m a(r). 

Such algorithm is reported and proved in [TB] ; later Mora in a series of lecture 
notes labelled the restriction of this decomposition and interpolation to the case 
of a set of distinct points as 'Axis-of-Evil' theorem and gave a precise descrip- 
tion, but no simple proof, of the result stated in [16J : S. Stcidel implemented 
the procedure in Singular [B], [18] . 

We give here such explicit algorithm that, given a finite set X of distinct points, 
provides a complete decomposition X = |Jm=i Ua=i ^>ms{T) on which, applying 
Cerlienco-Mureddu algorithm and interpolation, produces the required linear 
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factorization for a lexicographical minimal Groebner basis F = {fi,...,f r } of 
the ideal /(X) and thus a very simple proof of the 'Axis-of-Evil' theorem in this 
particular situation. 

This algorithm arranges the r terms ti belonging to G(J(X)) with respect to lex 
(ti < ••• < t r ) and constructs the factorization of each fi £ F through a suit- 
able interpolation on a subset S m $(ti) of X depending on the exponents of the 
corresponding t{. More precisely, Cerlienco-Mureddu give an algorithm that 
enables to find the Groebner escalier N(I(X.)) and the minimal basis G(/(X)) 
of the monomial ideal T(7(X)). 

Since the 'Axis of Evil' algorithm's starting point are the elements of X and 
the monomials of the finite Groebner escalier N (computed using Cerlienco- 
Mureddu algorithm), but the algorithm requires as input the monomial basis 
of T(/(X)), we also define the 'potential expansion' 's algorithm. 
It takes N and computes the minimal basis. 

1 note here that Marinari-Mora explicitly deduced, as trivial corollaries of their 
'Axis-of-Evil' procedure, Lazard theorem ([8]), Elimination theorem ([2]), Kalk- 
brener theorem ([13]), part of Gianni-Kalkbrener theorem ([7],[T2]); they how- 
ever remarked that, having being strongly influenced by Gianni-Kalkbrner re- 
sult, they cannot dismiss the possibility that Gianni-Kalkbrenner argument is 
an essential tool of their proof of the 'Axis-of-Evil' theorem. 

2 Notation. 

Let P := k[x%, ...,x n ] = (§) de ^Pd be the ring of polynomials in n variables and 
coefficients in the base field k. For all M C P, Md = M U Pd is its degree d 
part. Call T the semigroup of terms, generated by the set {xi, ...,x n }'- 

r-K 1 •••<», K-,0 en 

Letting a = (a%, a n ) € N", we will often write x a instead of x" 1 ■ ■ ■ x" n . 
Define also the set 

T[m] :=Tnk[x u ...,x m ] = {xf ■ ■ ■ x^ / {a u ...,a m ) G N m }. 

For each semigroup ordering < on T (i.e. a total ordering such that t\ < 
ti => tt\ < U2, Vt,t\,t2 € T) we can represent a polynomial / € P as a linear 
combination (with coefficients in k) of monomials arranged w.r.t. <: 

s 

f = J2 c (-f' < ) < = J2 : **) e fc *' u e r ' > - > ts - 

teT i=i 

We will call T(f) = Lt(f) := t x the leading term of / and tail(f) = f - T(f) 
the tail of /. 

We can also express it in a unique way as 

/ = ^g t x l n G k[xi, ...,x n -i][x n ], g t e k[x%, x„_i], g s ^ 

i=0 

(where S := deg n (f) is the degree w.r.t. x n ). 

We denote Lp(f) := gs, the leading polynomial of /. 
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Definition 2.1. For each monomial t £ T and Xj\t, the only u e T such 
that t = XjU is called j-th predecessor of t. 

A subset N C T is an order ideal if 

(el^se NVs\t. 

Let JV C Tan order ideal, A subset N C T is an order ideal if and only if 
T \ N = J is a semigroup ideal (i.e. t £ J tr £ J, Vt £ T). 
We set iV(J) :=T\T(J) = N. 

For a semigroup ideal J, G(J) denotes its minimal basis and 

G(J) := {r e J | each predecessor of r e iV( J)} = 
= {r e T | JV(J) U {r} order ideal, r £ iV(J)}. 

For all subsets G C P, we define T{G} := {T(g), g £ G} and we call T(G) the 
semigroup ideal {rT(g) 7 r £T,g £ G}, generated by TIG}. 
For any ideal I<P consider the semigroup ideal T(I) = T{I}, denoting by abuse 
of notation G(I) its minimal basis G(I) and the border ideal of / 

B{I) := {x h t, l<h<n,t£ N(I)} \ N(I) = 
= T{I) n ({1} U {x h t, l<h<n,t£ N(I)}). 

We will always consider the lexicographic order induced by x\ < ... < x n , i.e: 

xl 1 ■■■ x a n n < x\ x ■ ■ ■ x b n n ^ 3j | a,j < bj, a,i = hi, Vi > j. 

This is a term order, that is a semigroup ordering such that 1 lower to every 
variable or, equivalently, it is a well ordering. 

Lemma / Definition 2.2. We have: 

1. P = I®k[N{I)}; 

2. P/I^k[N(I)}; 

3. V/ £ P, 3\g := Can{f,I) = Eteivm 7(/,*, <)* e k[N(I)], called canoni- 
cal form of f with respect to I, such that f — g £ I. 

Definition 2.3. Given a term order < on T: 

1. a Groebner basis of I is a set G C I such thatT(G) = T{I}, that is 
T{G} generates the semigroup ideal T{I) = T{I}; 

2. a minimal Groebner basis is a Groebner basis such that divisibility 
relations among the leading monomials of its members do not exist; 

3. the unique reduced Groebner basis of I is the set: 

Q{I) := {t - Can( T , I) : r £ G(I)}. 

Each member of the reduced Groebner basis has a leading term which 
does not divide any monomial of another member. 
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Let X = {Pi, Pjy} C k n be a finite set of distinct points 
Pi := (aa,,..,a in ), i = 1,...,N. 

We call 

/(X) := {/ G P : f(Pi) = 0, V*}, 

the ideal of points of X. 

Finally, we define the projection maps: 

7r m : k n -)■ fc m 7r m : fc n -> fc"- m+1 

(Xi, .., X n ) i ^ (Xi, X m ), (Xi, .., X„) i V (x m , X n ) 

and, for P € fc™, X C fc", let 

n s (P,X) := {Pi G X/7r.(P,) = tt s (P)}, 

n s (p,x) := {p g x/tt s (p) = ^ S (P)}, 

extending in the obvious way the meanings of 7r s (d), 7r s (d), n s (d, D), IF (d, D) 

to d G N" c k n e P c N n C N™. 

With the same notation 7r m we indicate also 

3 Cerlienco-Mureddu Correspondence. 

Cerlienco and Mureddu ([3], [J], [5]) provided an algorithm which solves the 
following 

Problem: Given finite ordered set of distinct points 

X := (Pi,...,Pv) C k n - P := { ail ,...,a in ) 

compute a monomial basis (w.r.t. the lexicographic order induced by x\ < ... < 
x n ) of the quotient k[x%,..., cc„]//(X), where X denotes the support {Pi, Pn} 
ofX. © 

More precisely, they 

• define the operator associating to X an ordered set 

*(X):= (d 1 ,...,d A r)cN' i 

such that |$(X)| = |X| = N and such that , for all m < N the subset 
(di, ..d m ) is exactly $((Pi, ...,P m )). 

• define the a-value w.r.t. X s = er(P, X) of a point P G K n \ X as 
the maximal integer such that II s _i(P, X) ^ (by convention, VP, X, 

n o (p,x)^0). 

For P ^ X, they define the set 

S(P,X) := {P G X/tt s _i(P) = 7r s _i(P), s = a(P,X)} 
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containing all the points of X having the first s — 1 coordinates equal to 
those of P £ X. They extend the notation to the case P = Pj 6 X in the 
following way: 

a(P,X) :=a(P,{P 1 ,..,F J _ 1 }) 
E(P,X) :=E(P,{Pi,..,.Rf-i}). 

Remark 3.1. Given a term order a monomial basis for A := k[x%, x n ]/J(X), 

[x !l ], [x iN ], with x 11 X ... X x iN 

is called minimal wi£A respect to the term order if, for every other monomial 
basis [x'i], [x'n], with x 1 * X ... X x 1 ™ /or t/ie ^4 if holds 

Vj = 1, ...,7V, x ij ^ x'j. 
/n |3]/ ; they state that the computed monomial basis is the minimal one. 

Proposition 3.2. (J3j/) 

Let D := $(X). Then {[x d ]/d G D} is a monomial basis for k[x\, x n ]/I(X). 
Such a monomial basis is minimal with respect to the given <. 

Once the Groebner escalier N is known, it is very simple to compute the minimal 
basis G of T(J(X)) = T\N. Given the set X, the first step to compute the linear 
factorization of a minimal Groebner basis will be to apply Cerlienco-Mureddu 
algorithm to X and compute N, in order to obtain G. 

4 The potential expansion's algorithm. 

Consider the polynomial ring k[x\, x n ] with usual ordering <. Given a finite 
set of distinct points X = {Pi, Pat}, consider the ideal /(X) < k[x±, x n ] 
which is radical and zerodimensional, so its Groebner escalier N is a finite set. 
We will compute the minimal monomial basis G of T(J(X)), given the Groeb- 
ner escalier. The algorithm actually provides correct results irrispective of the 
given term ordering, but since we use Cerlienco-Mureddu correspondence, we 
will consider only our lex order. 

In order to make the reasoning clear, we will represent the monomials using the 

same diagrams introduced in |15) to study properties of Borel ideals. 

Apply Cerlienco-Mureddu correspondence to X in order to have iV(X) = 

{ti, ...,Tjv}. It is well known (see, for instance [16]) that |JV(X)| = |X|. 

We first define the potential expansion of a subset H C T, from which the 

algortihm bears its name. 

Definition 4.1. Let H C Tj for some j € N* we set C^(H) := H and, for all 
I G N* CW(t) = T j+ i \ {x u ...,x n } ■ (75 + i_i \ CV-V(H)). 

We then slice the Groebner escalier by degree, having No, N\, ■ ■ ■ Nh, where 
h is the maximal degree of terms appearing in TV. 

The minimal monomial basis G(/(X)) will have at most degree h + 1. As a 
matter of fact, if r G G with deg(r) = d > h + 1 its predecessors will belong to 
N and have degree d — 1 > h + 1 which is impossible. 
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Algorithm 1 Ccrlicnco-Murcddu algorithm. 
1: procedure CeMu(X) -> $(X) 
2: if N = 1 then 
3: *(X):={(0,...,0)}. 
4: end if 

5: if 1< N theni> suppose to know by induction hypothesis &((P± , . .. , P N _ i )) = (d^ , . . . , djy _ i ) 

and look for d N = &(P N ). 

6: s = a(P N ,X). 

7: for i = n to 1 do 

8: if i > s then 

9: d Ni = 0. 

10: end if 

11: if i = s then 

12: m, (1 < m < n), maximal s.t 7r s _i(P TO ) = 7r s _i(P/v), 

7r S ^ 1 (d m ) = (0, ...,0) = 7T S ^ 1 (d7v)- > Pm is the a -anteceiicnf of Pjy w.r.t. (Pi , . . . , P N _ j ) , 
*((P 1 ,...,P JV _ 1 )). 

13: dj\[ s — d ms + 1. 

14: end if 

15: if i < s then > wo use induction here. 

16: W(Pjv,X) := {P G X|$(P) = d = (*,... *,d JV8 ,0,...,0), } - 

{Pjl, Pjr}- 

17: Q:=n s ^(W(P N ,Xj). > 

Q| = W(P N ,X)I = r<N. If h < r = |W(P N ,X)|, then ir s _ 1 (P jft ) ^ 7r 3 _ 1 (P JV ). Moreover, since * is 
inductive, if h < k < r then ir 3 _ 1 (P j7l ) ? tt s _ i ( Pj fc ) . 
18: 7T s _l(d7v) =d r . > By the induction hypothesis,*^) = (di, ..,d^) and 

VI < i < r, 5; = TTg.itdji). 

19: break. 
20: end if 

21: end for 

22: end if 

23: return <f>(X). 

24: end procedure 
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The computation of G is performed as follows. 

Consider % Vi = 0, h+1: it is well known that \%\ = ("+17 1 ). 

For each i, define Gerii(I) := {t G G(I) \ deg{t) < i}. Since 7 is a proper ideal, 

Geno(J)=0. 

Let ft- the minimal i such that Gerih(T) 7^ 0, Vi > 1 

h+i-l 

Gen i+h = Gen h+i - 1 U(T h+i \(N h+i U (J C^ h+i ~ j \Gj))). 

j=h+i 

We then have 



Algorithm 2 The potential expansion's algorithm. 



1: procedure PotExp(7V(7)) -> 7 £> 7 is expressed using its minimal basis. 

Require: TV = [7V , N h , N h+1 ], such that iV^+i = 0. 

2: C = [0] • ^ the potential expansion' s list. 

Gen = 0. 
7 = (0). 

for z = to h + 1 do 

d=( n+de9 ^ ) - 1 )-\N i UC[t\\. 
if of = then 

C[z + 1] = PotentialExpansion(C[i]) . 
Gen,, = (0) 



t> no new generators. 



else 

Gem =Ti\{NiU C\i}). 

C[i + 1] = Potential Expansion(Gerii U C[i\). 
7 = 7 + Geni. 
end if 
end for 
return 7 
end procedure 



D> adding new generators. 



The algorithm uses a subroutine Potential Expansion such that 
PotentialExpansion(F) = C (1) {F). 

We will also have a subroutine finding T h+l \ {N h+i U Uj=fc+i C^ h+i -^{Gj)). 

WLOG we will think that the sets 7; +l and N h+t U Ujih+i C^+^iGj) are 
ordered w.r.t. the same ordering, since it is enough to perform a merging with 
the Groebner escalier and the potential expansion previously ordered. 
All these steps end: the subroutine finding the complementary can be developed 
performing a loop on the two ordered lists A := % = [ai,..,a m ], m > n and 
B := Ni U CW = [61, 6„] (using two indices i, j), keeping in mind that B <Z A 
m B = A and that 7? [7] > A[i] at every step. Start with b\\ if &i = a\ we set 
i = i + 1; j = ? + 1. 

If we find dj 7^ &j for a certain couple we put A[i] in the complementary 

and i = i + 1 without modifying j. 

Example 4.2. T/iere are situations in which N contains monomials of degree 
at most h, but also the minimal basis shares the same property. 
Take I = (x 3 , y 2 ,z 2 , xy) < k[x, y, z), whose Groebner escalier is: 
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A^o = {1} 
Ni = {x,y,z} 
N 2 = {yz,xz,x 2 } 
N 3 = {x 2 z}: 




1 

The monomial basis does not contain elements of degree 4. 

We call Gi the set of i-degree elements of the minimal basis and I the mono- 
mial ideal we want to find. 

Lemma 4.3. For all i = 0, h+1 

i-l 

T i \(N i u\JC i - j (G j )) = G i . 

3 = 1 

Proof: The inclusion T l \ (Ni U [fr} t C l -\G )) D d is trivial, so we only 

prove Tj \ (JV, U U}=i ^(Gj)) Q G l . 

Consider r G %\ (Ni U U^l C^iGj)). Clearly t e I. 

Let a the ith predecessor of r; if a e /, 3t e G with er = i ■ m for a suitable 
m e T- 

Then r = t • m ■ x t i.e. r e (Jj-=i ^(Gj). o 

This lemma assures that the result obtained via the potential expansion's algo- 
rithm is correct. 

5 The Axis of Evil Algorithm. 

A 0— dimensional radical ideal / < P is completely determined if we know the 
set V(I) of its zeros. 

Consider a finite set of distinct points X = {Pi, ...,P r }; we will denote indiffer- 
ently the Groebner escalier of the ideal /(X) with N(I(X.)) or N. A variation 
of Cerlienco-Mureddu algorithm ([3J) allows us to find a 'linear factorization' 
for every element of a lexicographic minimal Groebner basis in the sense of the 

Theorem 5.1. Let ti := xf 1 ■ ■ ■ x„" , i = 1, r be the generators of the minimal 
basis ofT(I), where I is a 0— dimensional radical ideal. 

A combinatorical algorithm and interpolation allow us to deduce polynomials 

'fmSi — X m graSi (x\ , . . . , X m — \ ) , 

V«, m, 5, with l<i<r,\<m<n,\<5< d m such that 

/* = n n imM ^ i 

m S 
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where /j, i = 1, ...,r are i/ie polynomials forming a minimal Groebner basis of 
I with respect to the lexicographic order induced by x\ < ... < x n . 

In that algorithm we will use the projections, as we defined in section [3J 
The Axis of Evil algorithm works then in the following way: 

• consider Tj :— xf 1 ■ ■ ■ x„ n € G. The required polynomial / = tj + tail(f) 
is factorized in Yl7=i ^» factors: d\ polynomials whose leading term is x\, 
d 2 polynomials such that their leading term is X2 and so on; 

• consider the monomials x^x 2 2 ■ ■ ■ x^ n such that ai < d%; 

• every such monomial is associated, via Cerlienco-Mureddu Correspondence, 
to a point of our set X. Project these points with respect to the first co- 
ordinate, obtaining d\ numbers y\, ...,y c i 1 ; 

• X\ — yi, i — 1, d\ are the first d± factors; 

• construct the subset D 2 o of X containing all the points in which the prod- 
uct (xi — yi) ■ ■ ■ (xi — does not vanish. If it is empty then stop and 
consider the next monomial in G; otherwise continue as follows; 

• find the set N 2 {tj) of all monomials in T[2] such that x^x^ 2 < x^x^ 2 ; 

• split the elements of N2 (Tj ) with respect to the exponents of X2 and con- 
struct, via Cerlienco-Mureddu correspondence, the set 

{<$>-\vxt- s xf ■ --xt)/ve T[l], vx d 2 2 - 5 € N 2 ( Tj )} 

• intersect the previous set with D20, project the resulting set of points 
(A^s (rj)) with respect to the first two coordinates and apply Cerlienco- 
Mureddu Correspondence, obtaining a set E 2 srj \ 

• interpolate over A 2 g(Tj), finding d 2 factors whose leading terms are all 
equal to x 2 . The monomials of E 2 s T . are the ones appearing in such 
factorization; 

• update the set of points in which the current polynomial does not vanish 
and stop if it is empty; 

• repeat these steps letting all the variables vary one by one; 

• repeat all the steps for all r, G G. 

Remark 5.2. Given Tj = x^ 1 ■ ■ ■ x d " € G, every variable Xi will appear only di 
times in the execution of the algorithm. 

Remark 5.3. The sets N m (Tj) := {lj € T[m], Tj > wi^ 1 • • • x^™ € N} (in 
particular for m = 1 we have N\(tj) := {x\/ i < di}) are constructed in order 
to determine in which points it is necessary to interpolate. 
Since for \i > Tj the Cerlienco-Mureddu correspondence provides a point P^i 
such that 3k € {1, ...,n} : i^k(P^) — ^k(P^'), in order to obtain polynomials 
vanishing on all the poinst of X it is not necessary to interpolate in the whole 
^ 1 (N) as it suffices to consider only those corresponding to /1 € N with (i < Tj. 
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Algorithm 3 The Axis of Evil algorithm. 

l: procedure AoE(X, G(I(X)) := {n, r r }) -> R > i? contains a factorizcd minimal 

Groebner basis of /. 

Require: the elements G(/(X)) arc in increasing order w.r.t the lexicographical 
order w.r.t. X\ < ... < x r . 



2: R = <D 

3: for i = 1 to r do 

4: JVi(Tj) := {x\l i<d 1 } = {L0e T[l], Tj > CJX^ 2 

5: Arfo) :={$- 1 (ii^-i*')/i«ij}cX. 

6: Bi(tj) := 7Ti(Ai(Tj)) C fc. 

7: 71t, := llaeS!^)^! -«)• 

8: for to = 2 to n do 

9: CrriTj :== E[l/=1 7vt^ • 

10: AnO := {Pi G X/ ( mT] {Pi) £ 0}. 

11: if |-D m o| = then 

12: R = [R, C, mTj \. 

13: break. 

14: end if 

15: N m { Tj ) := { W G T[m], rj > w^ 1 ■ • • x^ € N}. 

16: for 5 = 1 to d m do 

17: A m5 (r,) := "■^)/« e T[to - 

l], vxf rT - 5 e ^v m (r,)} n An(5-i)(Tj). 

18: E mS (Tj) := ^{■K m {A mS {Tj))). 
19: 

such that 7 m 5 Tj (P) = 0, VP € A m5 (Tj). 

20: £ mS := ll v=1 1ut 3 [l d =l Jmdr- 

21: £>m«(Tj) := {P € X/£ ro5 (P) /0}CX 

22: if |D m j(rj)| = then 

23: fl= [_R,£ mi5 ]. 

24: break. 

25: end if 

26: end for 

27: 7„ lTj := Il«7m«T r 

28: end for 

29: end for 
30: return R. 
31: end procedure 
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Remark 5.4. The terms smaller than tj mentioned before are found releasing 
all the variables one by one. 

Imagine the monomials in k[x\, ...x n \ as points in k n , identifying every term to 
the n-uple of its exponents. So we can 'draw' them in a n-dimcnsional space and 
we can think our realeasings as an increment by one of the 'directions' where 
we can move there. 

We point out that N m (Tj) C Nh{rj) for m < h. 

If a; € N m (Tj), Tj > wx^ 1 ■ ■ ■ xfr e N ; as uxf^ ■■■xfr \x%£ ■■■xfr we have 
ux d h+i ■ • • xfr € N and 

, l!T dh+1 . . . r d " < T dm+1 . . . r d " < T 

At each step we will count out all the points in which the polynomial already 
vanishes and we will stop the computation when the current factorized polyno- 
mial vanishes on the whole X. 
We will see an example of it later. 

Remark 5.5. If the number of released variables is > 1, we also must split the 
obtained monomials regarding the exponent of the maximal variable. 
Consider then the loop on S and, in particular, the set: 

C m5 (T 3 ) := {<&- V^-^+i 1 ' ■ -< n )h e T[m - 1], vxt~ 5 € N m ( Tj )}. 

We intersect C m s(Tj) with the subset of X containing the points not vanishing 
the current factorized polynomial. 

We can easily notice that, performing the algorithm, we only compute the sets 
C m i ( Tj ) , • • • , C m d m {Tj ) , but in N m (Tj ) there are also monomials co = x" 1 • • • x°^~l x\ 
such that Tj > ujx'l^l ■ ■ ■ x^™ € N, which would be generated considering (5 = 0. 
They are not considered in the algorithm because they are related to monomials 
examined in the previous step: = ij 1 • • -x^Zl € iV TO _i, so the corresponding 
points have already been treated. Taking (5 = 0, ..,d m , the sets C m s(Tj) form a 
partition of N m (Tj) basing on the degree of x m . As a matter of fact, in order to 
have lj € N m (Tj) we must have tj > wif^ 1 • • • x^™, where wx^Y • • • x^" e N, 
then the exponent of x m will be the first checked in the lexicographic test and 
so it will be limited by d m . 

According to the values of this exponent, the ones associated to smaller variables 
will vary. 

Remark 5.6. At the beginning of the algorithm, we imposed the monomials 

j = l...,r to be in increasing order with respect <. The steps made by the 
algorithm on each Tj are totally independent both on those made and on those 
to be made on a monomial (it is indifferent whether j ^ k) belonging to G, 
so we will obtain the same factorizations even if we launch the computation on 
a list of unordered monomials. 

Clearly, the result of our computation is not the reduced Groebner basis of the 
given ideal, it is only one of the minimal Groebner bases but we can obtain 
the reduced Groebner basis via simple reduction. 

We decided to put the monomials in such an order because we want every 
polynomial to be reduced with respect to the 'previous ' ones. 
If fj is one of our resulting polynomials and Lt(fj) — Tj, the polynomials 
utilizable to reduce fj (the previous ones) must be necessarily all and only the 
ones having as leading terms elements in G lower than the given Tj . 
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The algorithm terminates because it works on: 

1. points in the finite set X; 

2. monomials r e G (they are in a finite number, [16]); 

3. a finite set of variables. 

Let us study the correctness of the algorithm. 

Lemma 5.7. The factorized polynomials obtained from our algorithm vanish 
on all the points of the set X. 

Proof: Suppose we want to construct j T with r = x" 1 • • • a;"". 

Let /i = x^ 1 ■ ■ ■ x^ n , corresponding to a point P p € X through Cerlienco- 

Mureddu Correspondence. 

Let fi < t, then at least one of the exponents of the variables appearing in [x is 
lower than the corresponding in r, say < ccj, so fi is linked to an element of 
iVj(r) and so it can, alternatively: 

• belong to A«(t) for some S; 

• be such that the corresponding point already annihilates the polynomial 
found. 

If fi > t (since r ^ N, it is surely impossible that r = /i) then there will be a 
point P^t such that 

7r J -(P M ) = 7r j (P^) ) 

corresponding to a fi/ < r. 

We then use // and we come back to the previous situation, o 



Corollary 5.8. The ideal generated by these polynomials is exactly /(X). 

Proof: By the previous lemma, the polynomials vanish on all the points of 
the set X and the equality comes out by reasons of multiplicity o 

The resulting polynomials form a minimal Groebner basis because: 

• they vanish on all the points of X; 

• their heads form exactly G(J(X)). 

Notice that we can obtain the current interpolating polynomial applying Moeller 
algorithm to the projection through n m of the points of the current A m s(r) 

(ED- 

Example 5.9. Let 

X:= {(4, 0,0), (2, 1,4), (2,4,0), (3,0,1), (2, 1,3), (1,3,4), (2,4,3), (2,4,2), (1,0,2)}. 



Pi 
Pi 
P 3 
Pi 
P5 



= (4, 0, 0) : it is a single point, so $({(4, 0, 0)}) = (0, 0, 0) 
= (2,1,4) :s = l,m=l, (1,0,0) 
= (2,4,0) : s = 2, m = 2, (0,1,0) 
= (3,0,1) : s = 1, m= 1, (2,0,0) 
= (2,1,3) : s = 3, m = 2, (0,0,1) 
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P 6 := (1,3,4) : s = l, m = 4, (3,0,0) 

P 7 := (2, 4, 3) : s = 3, m = 3, W = {(2, 1, 3), (2, 4, 3)}, t 7 = (0, 1, 1) 
P 8 := (2,4,2) : s = 3, m = 7, (0,0,2) 

P 9 := (1,0,2) : a = 2, m = 6, W = {(2, 4, 0), (1, 0, 2)}, t 9 = (1,1,0). 

Then N :— {1, X\, x%, x\, 23, x\ , X2X3, irf, x x x 2 } and so we can easily obtain 

G = {xf, x\x 2 , x\, X1X3, X2X3, X3}. 

T/ie monomials belonging to G are exactly the input for the Axis of Evil algo- 
rithm and they are already ordered with respect to our ordering: 
starting with r = x\ we obtain 

N 1 {t) = {1,x 1 ,x\,x\}; 
A 1 (r) = {(4,0,0),(2,l,4), 
(3, 0, 1), (1, 3, 4)} : these are the corre- 
sponding points via Cerlienco-Mureddu 
Correspondence; 
B x {r)= {4,2,3,1} 

7 i T = {x x - 4)(xi - 2){xi ~ 3)(n - 1); 

i/ie linear factors are only depending 
from x\ are computed in the same time, 
m = 2: 
(2t = 7Xr 

-D2o( T ) — 0.' stop here obtaining, as first result, a polynomial having as lead- 
ing term an element of G (while the other monomials belong to N ) and belonging 
to /(X) since it vanishes in every point of X (so an element of our minimal 
Groebner basis). 
t = x\x 2 Ni(t) = {l,xi}; 
A 1 (t)={(2,4,0),(1,0,2)}; 
Bi(r)={2,l} 
Tlr = (xi - 2)(xi - 1) 

m = 2: 

Am,(t) = {(4,0,0),(3,0,1)} 
We can not stop here, since the ob- 
tained polynomial does not vanish at all 
the points and its head is different from 
t € G. 

* * N 2 (t) = {l,xi,x\,x\,x 2 ,xix 2 }; doing 

so, we find all the monomials of the 
previous step (we computed their corre- 
sponding points) and some new ones. 
6=1: 

A a i(r) = {(4,0,0) J (3,0 J l)} = U 20 
The monomials vx^ 5 are 1 l , corresponding to the points Pi , P 2 , P4 , Pq . 

The polynomial already vanishes on P 2 ,Pg, so we consider only the remaining 
two points. 

E 21 {t) = {l,xi}. 

721r = x 2 ; 

61 = 7it721t = (xi - 2)(.ti - l)x 2 ; 
D 21 {t) = 0. 





13 



Remark that ^ 2t is actually j 2 i T . 

2 

T = x 2 

N 1 {t) = %; 
A x (t) = 0; 
Bi(r)=0 

m = 2: 

#20 (r) = X 



-W 2 (t) = {l,a;i,xf,xf,X2,a;ia;2}; 5 = 1: 

A 21 (t)={(2,4,0),(1,0,2)}; 

E 21 {t)={1,x 1 }; 

721r = x 2 - 4xi + 4 

6i = 7ir72ir = «2 - 4xi + 4; 

£> 21 (t) = {(4,0,0), (2, 1,4), (3,0,1), 

(2. 1.3) , (1, 3,4)}; 
6 = 2: 

Aaa(T)={(4,0,0),(2,l,4),(3,0,l), 

(1.3.4) } 

The terms raJ^T 5 are l,x\,x\,x\ 
and they correspond exactly to 
Pi, P2, Pi, 



E 2 2(r) = {l,a;i,a;f,^}; 

722r = 2x 2 — arf + 7xi - 12; 

62 = (x 2 - 4xi + 4) (2x 2 -xj + 7xi - 12) 

D 22 (r) = 0; 

T = X1X3 



Ni(r)={l}; 
A 1 (r) = {(2,l,3)}; 
Bi(r) = {2} 

7lr = (#1 - 2) 



m = 2: 
^2(r) = {1}. 

D 20 (t) = {(4, 0, 0), (3, 0, 1), (1, 3, 4), (1, 0, 2)} 
8=1: 

D 21 (t) = D 2Q (t); 
m = 3: 

^3(7") = {l,xi,x 2 ,Xi,X3,x?,xix 2 }; 

Cmr = (Xl - 2); 

^3o(r) = {(4, 0, 0), (3, 0, 1), (1, 3, 4), (1, 0, 2)}; 



E 3 i{t) = {l,xi,x^,x 2 }; 



6=1: 

A 31 (t) ={(4, 0,0), (3,0,1), (1,3,4), 
(1,0,2)} 

The terms are 1, xi, x\ , x\ , X2, X1X2, 
corresponding to P\, P 2 , P3, Pa, Pq, 
Pq , but we can neglect P 2 ,P 3 . 
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731 (t) = 6^3 - 4a; 2 + x\ - x\ - 12; 
£ 3 i = (xi - 2)(6x 3 - 4x 2 + x\ - x x - 12); 
Att(r) = 0. 
The desired polynomial is 73 T = 731 (t). 

T = X 2 x\ 

N 1 (t)=%: 
Ai(t)=0; 
Bi(t)=0 
m = 2: 



iV 2 (r) = {l}; 
D 20 (t) = X; 
<J= 1: 

A 21 (t) ={(2,4,2)}; 
£ 2 i(r)={l}; 

721r = 2>2 - 4 

61 = x 2 - 4; 

£» 21 (r) = {(4,0,0),(2,l,4),(3,0,l), 
(2, 1,3), (1,3, 4), (1,0, 2)}; 



m = 3: 

C3r = X 2 - 4 

Aso(t) = £> 2 i(r); 



N 3 (t) = JV(X); 
(5 = 1; 

A 31 (r)= {(2,1,3)}. 
£ 3 i(t)={1}; 

721r =3^3-3 

&i = (a:2-4)(a;3-3); 

£» 3 i(r) = {(4,0,0),(2,l,4),(3,0,l), 

(1,3,4), (1,0,2)}; 



5 = 2: 

A 32 (t)=D 31 (t); 

E 32 {t) = {l,x 1 ,xj,x 3 1 ,x 2 }; 

732r =x 3 - 4x 2 - + Alxj - 96xi + 48; 

62 = (s a - 4) (a* - 3) (a* - 4^2 - 5x? + 41a:? - 96xi + 48); 

£32 (r) = 0; 

73r = (2:3 - 3)(x 3 - 4x 2 - 5x? + 4lxj - 96xi + 48); 

T — x l 
#i(r) = 0; 
Ai(t)=9; 
Bi(t) =0 
m = 2: 
D 20 (t) = X; 
N 2 (t) =0; 
5=1: 
A 2 i(r) = 0; 
£>2i(r) =X; 



15 



m = 3: 
#30 = X; 

I 



N 3 (t) = JV(X); 5 = 1; 
^(r) ={(2,4,2)}; 
#3i(r)={l}; 

• > 73ir = £3 - 2; 

• 61 = £3 - 2; 

£> 3 i(r) = {(4, 0,0), (2, 1,4), (2, 4,0), 
(3,0,1),(2,1,3),(1,3,4),(2,4,3)}; 
6 = 2: 

A 32 (t) = {(2, 1,3), (2, 4, 3)}; 

£ 32 (r) ={W; 

732r = - 3; 

62 = (as 3 - 2) (a; 3 - 3); 

# 32 = {(4, 0, 0), (2, 1, 4), (2, 4, 0), (3, 0, 1), (1, 3, 4)}; 
5 = 3: 

^33 (t) = -D32; 

E 33 (t) = {l,x 1 ,xl,xl,x 2 }; 

733r = 6x 3 + 8x 2 - 5x? + 35xj - 54xi + 24; 

£ 33 = (13 - 2)Os - 3) (6a; 3 + 8x 2 - 5xf + 35x1 - 54xi + 24); 

D 33 (t) = 0; 

The required polynomial is 7 3r = (x 3 — 2) (213 — 3)(6x 3 + 8x2 — 5x 3 + 35x\ — 
54a;i+24). 

Then our minimal Groebner basis of the ideal associated to X with respect to 
the given order is: 

Q(I(X)) = jxf - IOie? + 35x1 - 50xi + 24, x 2 x\ - 3x 2 xi + 2x 2 , 

x\ — 2x2X1 — X2 + 2x 3 — 16xf + 38xi — 24, x 3 x — 2x 3 — |x2Xi + |x2 + 

+^x 3 - \x\ - fxi + 4, x\x 2 - 4x§ - 7x 3 x 2 + 28x 3 + |x 2 xi + 
+f x 2 - fx 3 + 48x 2 - S^xi + 32, x 3 - 5x| + |x 3 x 2 - ^x 3 - f x 2 x x 

_40„ , 73^,3 197 2 i 1358 79 1 

— 9^X2 -t- — x 1 ~ x i 1 ~~ ' z p 

obtained by our polynomials by the reductions stated in the Axis of Evil Theorem. 
Finally, we remark that: 

1. let Tj = xf 1 ■ ■ ■ xfj" £ G. The polynomial we are looking for must contain 
exactly Y2?=i factors. It is impossible that the algorithm stops before, 
so it is impossible that a partial product vanishes on the whole X. In fact, 
if so, there would be a polynomial f £ I such that T(f) £ (G) (we know 
the minimal basis G before starting the Axis of Evil process) ; 

2. if we otain a factorized polynomial / such that its leading term T(/) 
belongs to the minimal basis G, then / vanishes over all X, because of 
15771 
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Example 5.10. Consider the following ideal, given with its primary decompo- 
sition: 

J := (xf,x 2 + X\, x 3 ) n (x\,x 2 - Xi,x 3 - 1) = 

= (xf,xix 2 ,xl,xix 3 - \xi - \x 2 ,x 2 x 3 - \xi - \x 2 ,x\ - x 3 ) <C[xi,x 2 ,x 3 ]. 

Call its generators fi,..., fe, considering them in the correct order. 

It is ^-dimensional because x\,x\,x\ <= In(J) (see fl&j), but it is not radical: 

its radical is \fj — (x 2 ,x\ — x 3 ,Xi). 

For such an ideal the Axis of Evil does not hold. 

Consider the polynomial fi = xix 3 — \x\ — \x 2 . 

By the Axis of Evil theorem \5.1\) . its factorization should be of the form: 

{xx + ...){x 3 + ...) 

and we should have 

xix 3 - -xi - -x 2 + Pxl + Qx\x 2 + Rx 2 , P,Q,Re C[xi,x 2 ,x 3 ], 

since we can only reduce deleting the multiples of x\, x\X 2 , x\, in order to obtain 
fi. In order to have the correct product we must have —\x 2 in it. We can not 
obtain it through reductions, so the only chance is that we have a product of the 
form 

k * hx 2 , 

with h,k constants such that hk — — |, in particular both different from 0. 
A priori, we can have two possibilities: 

• (xi + k + ...)(x 3 + hx 2 + ...); 

• (xi + hx 2 + ...)(x 3 + k + ...). 

The second one is impossible: the polynomial having x\ as head can not contain 
variables greater than x\, so we consider only: 

(xi + k + ...)(x 3 + hx 2 + ...). 

We will then obtain 

xix 3 + hxix 2 + kx 3 - ^x 2 + ... 
We can delete the term X\X 2 but it remains kx 3 which can not be reduced. 

6 Corollaries. 

We enumerate here some famous theorems which can be easily proved as corol- 
laries of the Axis of Evil Theorem. For more details see, for example, [IB] . 
Here we provide the general statements of these results, but clearly they can 
only be deduced under the hypothesis of the Axis of Evil theorem 
The first one is hazard Structural Theorem, which describes the structure of a 
minimal lexicographical Groebner basis of an / < k[x\, x 2 \. 
The original proof considers P = k[x\,x 2 ] = k[xi\[x 2 ] and it is based on the 
fact that k[xi] is a Principal Ideal Domain (PID). 

Norton-Salagean [T7] reformulated it using, more generally, R[x] with R PIR. 
We briefly recall the following 
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Definition 6.1. The content r/ £ R, with R PIR, of a polynomial f(x) £ R[x] 
is the GCD of its coefficients. A polynomial f{x) £ R[x] is called primitive if 
r f = 1. 

The primitive part of f{x) £ R[x] is the polynomial po(x) £ R[x] such that 

f(x) = r f p (x). 

Let R be a PIR, P := R[x]. Let I<P e F := {/q, f s } a minimal Groebner 
basis of I ordered in such a way that, called d(i) := deg(fi), Vi, < i < s 

d(0) < ... < d(s). 

Define then a — lc(fi), r, € i?\ {0} e pi £ P the leading coefficient, the content 
and the primitive part of fi, for all 1 < i < n. 

Theorem 6.2 (Lazard). I/ 7 moreover, R is a PID, then: 

• Jo = fGj • • ■ G s +i; 

. f j = PH j G j+1 >>-G a+1 , 1< j < s. 

L d(l) < ... < d(s); 

I?. Gi € i?, 1 < i < s + 1 is suc/i £/ia£ Cj_i = GjCj 

3. P = po (the primitive part of fo £ i?[x] 

4- Hi £ i?[x] is a monic polynomial of degree d(i) in x, for all i; 

5. for all i we have Hi + i £ (G\ ■ ■ ■ Gi, H\Gi ■ ■ ■ Gi, Hi^iGi, Hi); 

6. n = G i+ i ■ G s 

Theorem 6.3 (Norton-Salagean). With the previous notation, each 

Pi € (fj,3 < i) ■ n- 

In fact, we have = l\"n=i I\s=i lrn.su and pi = n«=i InSu- 
The second well-known result which can be straightforwardly derived from the 
Axis of Evil Theorem is the well known Elimination Theorem (see [5] for details) 

Theorem 6.4 ([15]). Let I <k[ ideal, take the lexicographical or- 

dering induced by x\ < ... < x n and call Ij the j-th elimination ideal Ij = 
I n k[x\, Xj\. Let Q be a Groebner basis ofL, then Qj = Q D k[xi, Xj] is 
a Groebner basis of Ij . 

The following result, Kalkbrener theorem (|13|. |16j). is another consequence 
of the Axis of Evil Theorem and it is a stronger characterization of the lexico- 
graphical ordering. 

For each subset L c k[x\, x n ], i — 1, — , n, VS £ N set 

Lis = {p £ L, \p £ k[xi, ...,Xi], degi(p) < 6} 

and 

Lp iiS = {Lp(p), p £ L i;S }. 
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Theorem 6.5 (Kalkbrcnner) . With the previous notations, considered an ideal 
J< k[xi, ...,x n ] and a Groebner basis Q of it, these forms are equivalent: 

• Q is a Groebner basis of I w.r.t, the lexicographical order < induced by 
X\ < ... < x n ; 

• Lpi t s(G) is a Groebner basis of Lpi t s(I), i = l,...,n, V<5 € N. 

Let us now mention Gianni- Kalkbrener theorem, whose situation is a bit 
more complicated (see [12], 0, [IS])- 

Theorem 6.6 (Gianni-Kalkbrcner). Let I < k[x\, x n \ an ideal and Q w.r.t 
the lexicographical order < induced by x\ < .... < x n . As before we define also 

g d = gn k[x x ,...,x d ]. 

Consider a = (pi, bj) G V(Id) and define the projection map 
$ Q : k[xi, x n ] k[x d+1 , x n ] 

f(xi, ...,x n ) i-> f(b\, ...,b d ,x d+1 , ...,x n ). 

Let a be the minimal value such that <$> a (Lp(g a )) ^ and j,S the values such 
that 

g a = Lp(g <y )x 5 j ^ 1 + ... € k[x!, ...,x 3 ] \ k[xi, ...,x 3 -i}. 

Then 

1. 3=5 + 1 

2. VgeGd, $ a (. 9 )=0; 

3. Vgeg d+ s, $a(.9) = 0; 

4- ^a(ga) = gcd($ a (g), g e G d +i) € k[x d+1 ]; 

5. V6 e k, (h, b 2 ,b) e V(I d+1 ) «■ $„(<?„)(&) = 0. 

Clearly (1 — 3) are essentially a corollary of theorem 16.31 on the other side, 
(4 — 5) apparently cannot be deduced from the Axis of Evil Theorem. 
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